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Abstract, imparity violating supersymmetric models (RPV SUSY) are becoming increasingly 
more appealing than its R -parity conserving counterpart in view of the hitherto non-observation 
of SUSY signals at the LHC. In this paper, RPV scenarios where neutrino masses are naturally 
generated are discussed, namely RPV through bilinear terms (bRPV) and the “/i from z/” 
supersymmetric standard model (iiv SSM). The latter is characterised by a rich Higgs sector that 
easily accommodates a 125-GeV Higgs boson. The phenomenology of such models at the LHC 
is reviewed, giving emphasis on final states with displaced objects, and relevant results obtained 
by LHC experiments are presented. The implications for dark matter for these theoretical 
proposals is also addressed. 


1. Introduction 

Supersymmetry (SUSY) [T is an extension of the Standard Model (SM) that assigns to each SM 
field a superpartner field with a spin differing by half a unit. SUSY provides elegant solutions 
to several open issues of the SM, such as the hierarchy problem, the nature of dark matter |2j, 
and the grand unification. SUSY is one of the most relevant scenarios of new physics probed at 
the CERN Large Hadron Collider |3j, yet no signs of SUSY have been observed so far. In view 
of these null results in conventional SUSY searches, it becomes mandatory to fully explore non¬ 
standard SUSY scenarios involving R -parity violation (RPV) [4j and/or quasi-stable particles. 

R parity is defined as R = ( — i) 3 (^- L )+ 25 / where B (L) is the baryon (lepton) number and 
S the spin, respectively. Hence all Standard Model particles have R = +1, whereas R = — 1 
for all SUSY particles. It is worth emphasising that the conservation of R parity is an ad-hoc 
assumption. The only strict limitation comes from the proton lifetime: non-conservation of both 
B and L leads to a rapid proton decay. R -parity conservation has serious consequences in SUSY 
phenomenology in colliders: the SUSY particles are produced in pairs and, most importantly, 
the lightest supersymmetric particle (LSP) is absolutely stable and weakly interacting, thus 
providing the characteristic high transverse missing momentum (E™ 1SS ) in SUSY events at 
colliders. Here we highlight two RPV models: the bilinear RPV (bRPV) and the [i- from-z/ 
supersymmetric standard model (/xz/SSM), which —as we shall see— both reproduce correctly 
the neutrino physics observations. 

The structure of this paper is as follows. Section[2]is dedicated to bRPV SUSY models: their 
connection with neutrino physics is explained in Section 2.1, its phenomenology at the LHC 


is provided in Section 2.2, while Section [2.3| presents the current constraints from the ATLAS 
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experiment. The fi- from-z/ supersymmetric standard model is discussed in Section [3] After a 
brief review of the theoretical motivation for [w SSM, descriptions of its possible signatures at 
the LHC based on multileptons and displaced vertices follow in Section |3T} Sections |3.2| and |3.3| 
focus on Higgs decays and unusual Z/W boson decays in the context of /wSSM, respectively. 
Aspects of RPV SUSY linked to dark matter are discussed in Section [4j The paper concludes 
with a summary and an outlook in Section [5j 


2. Bilinear R- parity breaking 

R -parity conservation (RPC) has several consequences such as the stability of the LSP, which 
is a weakly interacting massive particle (WIMP) and consequently a candidate for dark matter 
(DM) |2[. Being a WIMP, once produced at the LHC, it will escape detection, resulting in 
large missing transverse momentum, E™ 1SS . Providing a DM candidate is one of the strongest 
arguments in favour of RPC SUSY, nonetheless RPV models do exist that can explain DM 
through, for instance, very light gravitinos [5-1 1 


axions 110.121 or axinos 110,13,14 . There is 


no fundamental reason for R parity to be conserved thus lepton and baryon number violating 
renormalisable terms may appear in the supersymmetric potential in the following way: 


W — XijkLiLjE^ + X'-frLiQjD^ + CiLiH u 


+ X " ijk UtD]D% 


(i) 


where the couplings A, A' and X" are 3x3 Yukawa matrices — i,j, k being flavour indexes—, ci 
are parameters with units of mass and Q, £/, D, L, E , H U: Hd refer to supermultiplets. The first 
three types of terms lead to lepton number violation, while the baryon number is violated by 
the fourth one. 

As long as the breaking of R parity is spontaneous, only bilinear terms arise in the effective 
theory below the RPV scale, thus rendering bilinear R -parity violation a theoretically attractive 
scenario. Moreover, the bilinear model provides a theoretically self-consistent scheme in the 
sense that trilinear RPV implies, by renormalisation group effects, that also bilinear RPV is 
present, but not conversely [TB] . 

In other words, the simplest way to break R parity is to add bilinear terms to the MSSM 
potential. Besides that, additional bilinear soft SUSY breaking terms are introduced, which 
include small vacuum expectation values for the sneutrinos. The relevant superpotential W and 
the soft supersymmetry breaking terms V^oft? which include bilinear R -parity violation, would 
then be 16 


W 


W MSSM + 


Voft = v 8 


MSSM 

soft 




( 2 ) 

( 3 ) 


where the have units of mass. In fact, if SUSY was not broken, the bilinear terms could 
be rotated away and be converted into trilinear terms, however the presence of the soft SUSY 
breaking terms BiCiLiH u gives bRPV a physical meaning 


16 


There are nine new parameters introduced in this model: e$, Bi and the latter being the 
sneutrino vacuum expectation values (VEVs). However, after electroweak symmetry breaking 
and taking into account constraints from neutrino experiment results, only one free parameter 
remains in the model, which is set to be of the same order as the others. 


2.1. Connection with neutrino physics 

Sneutrino VEVs introduce a mixing between neutrinos and neutralinos, leading to a see-saw 
mechanism that gives mass to one neutrino mass scale at tree level. The second neutrino mass 
scale is induced by loop effects 117,18]. The same VEVs are also involved in the decay of the 
LSP, which in this case is the lightest neutralino. This implies a relation between neutrino 
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physics and some features of the LSP modes. An example of such a connection is given by the 
relation 1171 


tan 0 a tm — 


2 BRjft -> 

A r BR(xi ->• r^T) 


BRjft -> 

BR(x i ->■ 'r ± W / ) ’ 


( 4 ) 


where 0 a t m is the atmospheric neutrino mixing angle and the “alignment” parameters A i are 
defined as A^ = fiVi + with the VEV of H d . This relation between RPV and neutrino 
physics allows setting bounds on bRPV parameters from results of neutrino experiment and 
astrophysical data 1191. In the opposite direction, a possible positive signal observed in colliders 
may lead to the determination of some of the bRPV phenomenological properties, which in turn 
can constrain neutrino-physics parameters 1201. 


2.2. Phenomenology at the LHC 

In the specific bilinear R -parity violating model discussed here, the LSP is the lightest neutralino, 
The bRPV terms are embedded in the minimal Supergravity (mSUGRA) model, which 
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Xi 

imposes some restrictions which reduce the large number of parameters of the MSSM. The 
number of parameters in mSUGRA is reduced to five, namely mo, the scalar mass; mp 2 , the 
gaugino mass; Ag, the trilinear scalar coupling; tan/3 = v u jv^ the ratio of the Higgs VEVs; and 
sgn /i, the sign of the higgsino mass parameter. 

The six bRPV parameters for each mSUGRA point are determined by the SPheno 1 22 1 
spectrum calculator. SPheno uses as input the mSUGRA parameters and the neutrino physics 
constraints, and delivers as output the bRPV parameters (together with the mass spectra and 
the decay modes) compatible with these constraints. Once those quantities are calculated for 
a given set of bRPV parameters within an mSUGRA benchmark point, they can subsequently 
passed to an event generator and produce collision events at the LHC. 

It is stressed that the sparticle spectrum for bRPV-mSUGRA is —within theoretical 
uncertainties— the same as in RPC mSUGRA; it is the LSP decay modes and its lifetime 
that depend on the bRPV parameters. Typical Xi decay modes are shown in Fig. [l] for a model 
where the parameters tan/3 = 30, /a > 0 and Ag = —2mg have been chosen so as to obtain a 
125-GeV Higgs boson in agreement with the ATLAS 23 and CMS 124] discovery. The neutralino 


decays are dominated by leptonic (e, p) or r channels, making lepton-based searches ideal for 
constraining this model. The decay to a muon and a hadronically-decaying W has been studied 
in detail and it has been shown that it may lead to the measurement of the Xi mass through the 

Furthermore the fact 
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reconstruction of the Xi TQ'q peak in the fijj invariant mass 
that in the low-ra^ high-mo region the Xi is slightly long lived, as evident from Fig. [ 2 J opens 


up the possibility to use searches for displaced vertices |27,28| in order to probe this model at 
Lastly, the Xi decays to one or two neutrinos give rise to a moderately high £?j? lss , 
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the LHC 

thus rendering some Elf lss -based analyses pertinent to bRPV, as we shall see in the next section. 

So far we have discussed the bRPV in the context of supergravity models. If the bRPV terms 
are introduced to anomaly-mediated SUSY breaking (AMSB) instead, the Xi LSP decay modes 
do not change, whereas the wino-like neutralino (as opposed to the bino-like in mSUGRA) is 
characterised by stronger interactions, hence it is easier to be produced at LHC. The mass 
degeneracy of the Xi an d the xt makes the X\ long lived, decaying dominantly through RPV 
couplings to Ml , rll , Ibb, rbb 30 . In other words, the common feature of bRPV phenomenology 
is the presence of delayed decays leading to interesting signatures at the LHC. 


2.3. Constraints from ATLAS 

As mentioned above, thanks to the abundant neutrino production from the LSP decay, bRPV 
events at the LHC are expected to feature relatively (when compared to Standard Model 
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Figure 1 . Xi decay modes for a 
bRPV-mSUGRA point. 


mSUGRA-bRPV: tanp = 30, ju >0, A q = -2m 0 



m o [GeV] 


Figure 2. Proper decay length cr in millimetres for the 
lightest neutralino LSP in bRPV-mSUGRA parameter 
plane (mo, ^ 1 / 2 ) an d tan/3 = 30, /1 > 0 and Aq = 
—2mo- 


processes) high missing transverse momentum. High lepton/r multiplicity is also expected 
from the LSP decays and from upstream lepton production in the SUSY cascade decay, if strong 
production is considered. Both features are exploited when looking for a signal of the bRPV- 
mSUGRA model, as demonstrated in Ref. [25], where a first detailed Monte-Carlo-based study 
on the observability of the Xi I^q'q a t the LHC has been carried out. The very first bounds set 
on a bilinear RPV model in colliders were provided by an inclusive search for high E,f lss , three 


or four jets and one muon at yfs — 7 TeV and 1 fb 1 of ATLAS data [26,31 . The Xi decay 


to muonic channels is enhanced when compared to electrons, as shown in Fig. [TJ hence muon 
channels yield stronger limits than their electron counterparts. A slight loss of sensitivity in the 
high-mo low-m!^ region was observed due to the large LSP lifetime, which causes signal muons 
to be rejected by the cosmic-muon veto, i.e. the cut on muons with high impact parameter. 
Aforesaid limits were extended further by an analysis with 5 fb -1 based on events with high jet 
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multiplicity (at least seven jets), large E^f lss and exactly one lepton (muon or electron) 

Further constraints on bRPV-SUGRA for parameters tuned to attain a mass of the lightest 
Higgs equal to 125 GeV have been set recently by ATLAS (34f 36[ with the full data set of 
~ 20 fb -1 recorded at y/s = 7 TeV. The LSP decays to taus have been used as a handle to 
constrain the model in an analysis combining event selections with at least one tau, two taus or 
one tau and a lepton (e or fi) in conjunction with large Elf lss 35 . The observed limits can be 


consulted in the left panel of Fig. [3] The bounds set by an analysis seeking two same-sign leptons 
or three 5-jets and jets ;34] are significantly extended with respect to the r-analysis. As shown in 
the right panel of Fig. |3| values of mp 2 are excluded between 200 GeV and 490 GeV at 95% CL 
for mo <2.1 TeV. Interestingly enough, in both analyses, signal points with mp 2 < 200 GeV 
were not considered because the lepton acceptance was significantly reduced due to the increased 
LSP lifetime in that region, as seen previously in Fig. [2j Slightly stringent limits have been set 


36 


by ATLAS very recently by requiring one hard lepton, jets and large E^f lss 

The fact that in a large part of the parameter space the Xi LSP decays to a Z boson 
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should 


and a neutrino implies that searches for events containing Z and large E^f lss 
be pertinent to explore bRPV as well. In fact this LSP decay —as much as decays to rW 
and 11 W — is considerable not only when bRPV is embedded in mSUGRA, but also when 
bRPV is introduced to other well-motivated SUSY scenarios such as “natural” SUSY (391 
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and pMSSM |401. Ongoing work continues in interpreting the results of the same-sign lepton 
analysis [34] in natural pMSSM 04] with bRPV |421. 


Bilinear RPV Model, A 0 =-2m 0 , tan|3=30, p>0 



m 0 [GeV] 


1400 


MSUGRA/CMSSM: tan((3)=30, A = -2m 0 , p > 0, bRPV 
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Figure 3. Expected and observed 95% CL exclusion limits in the bilinear R -parity violating 
model obtained by a r + E™ 1SS analysis (left) | 35 || and by the two-same-sign-lepton analysis 
(right) [ 34 ] . The latter results are obtained by combining the electron and muon channels. The 
band around the median expected limit shows the ±1 a variations on the median expected limit, 
including all uncertainties except theoretical uncertainties on the signal. The dotted lines around 
the observed limit indicate the sensitivity to ±1 a variations on these theoretical uncertainties. 


3. The /i-irom-v supersymmetric standard model (/xz/SSM) 


The /w SSM 43,44 is a supersymmetric standard model that solves the (i problem 45 of 


the minimal supersymmetric standard model (MSSM) using the R -parity breaking couplings 
between the right-handed neutrino superfields and the Higgs bosons in the superpotential, 
XiV^HfiH u . The /i term is generated spontaneously through sneutrino vacuum expectation 
values, /i — A i(v£) : once the electroweak symmetry is broken, without introducing an extra 
singlet superfield as in the case of the next-to-MSSM (NMSSM) [46] . The complete /ivSSM 
superpotential is given by 


w = e ab (Y Uij H b u Qi u] + Y dij Q\ d) + Y eij H% L\ e] 
+ Y Vij Hi - Xi Of H a d H b u ) + 


(5) 


The couplings f° r bid a global U (1) symmetry avoiding the existence of a Goldstone 

boson, and also contribute to spontaneously generate Majorana masses for neutrinos at the 
electroweak scale. The latter feature is unlike the bilinear RPV model, where, as mentioned in 
Section |2.1[ only one mass is generated at the tree level and loop corrections are necessary to 
generate at least a second mass and a neutrino mixing matrix compatible with experiments. In 
the bilinear model, the //-like problem |47] is augmented by the presence of three bilinear terms. 

The /ivSSM phenomenology is largely defined by the parameters 


A, v c , tan/3, Mi, A\, A K , 


( 6 ) 
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where A = v^3A it the singlet-doublet mixing parameter (if universal A i are assumed), k is the 
common A\,A K are the soft SUSY-breaking parameters and M\ is the U{1) mass scale. 

In the fi z/SSM, as a consequence of R -parity violation, all neutral fermions mix together 
into ten neutralinos, and five charginos, xt- Since the three lightest neutralinos are the 
left-handed neutrinos, the “true” LSP would be the x®. Likewise the three lightest charginos 
xf coincide with the three charged leptons. Similarly, all scalars mix into eight CP-even, 
and seven CP-odd neutral Higgs bosons, P^, mass eigenstates. The three lighter neutral 
scalars, = 1,2,3, are in fact naturally light singlet-like states, leaving the fourth one, S 4 , 
to play the role of the discovered SM-Higgs-like scalar. Charged scalars, on the other hand, 
form seven mass eigenstates, P±- Analyses of the p^SSM, with attention to the neutrino and 
LHC phenomenology have been addressed in Refs. [44,48 . Other analyses concerning cosmology 


such as gravitino dark matter and electroweak baryogenesis can be found in Refs. [49-511 and 
in Ref. |52], respectively. In conclusion, /ivSSM is a well-motivated SUSY model with enriched 
phenomenology and notable signatures, which definitely deserve rigorous analyses by the LHC 
experiments. Its enlarged Higgs sector can easily accommodate the observed 125 GeV Higgs 
boson 


23,241 


3.1. Signatures at LHC 

Here a collider analysis together with detector simulation of an intriguing signal in the /ivSSM 
featuring non-prompt multileptons at the LHC, arising from the beyond SM decay of a 125 GeV 
scalar into a pair of lightest neutralinos, x% is presented 

decays into a scalar/pseudoscalar ( h/P ) and a neutrino, with the h/P further driven to decay 
into t + t _ , giving rise to a 4r final state, as shown in Fig. [4| The small P-parity breaking 
coupling of X 4 renders it long lived, yet it decays inside the inner tracker, thereby yielding 
clean detectable signatures: (i) high lepton multiplicity; and (ii) charged tracks originating from 
displaced vertices (DVs). 


1531. Since R parity is broken, each x\ 



Figure 4. The fiis SSM decay chain studied 
in Ref. [53]. The Higgs boson is produced 
through gluon fusion and has a mass of 125 GeV. 
The neutralino X 4 is long-lived and gives rise to 
displaced r leptons. 


The fiv SSM is characterised by the production of several high-px leptons 53 . Electrons and 


muons are produced through the leptonic r decays, yet muon pairs can appear directly through 
hi/Pi decays as well. With the chosen decay mode, the r multiplicity is considerable even though 
the r-identification efficiency is much lower 50%) when compared to that of electrons and 
muons (> 95%). Occasionally highly collimated QCD jets can fake hadronic r leptons, Th a d> 
and, as a result, r multiplicity may exceed the expected number of four. This faking however 
disappears with a higher px cut, which should also be sufficient to provide a single-lepton trigger 
for the analysis. 

Apart from the requirement of at least three or four leptons (including taus), a high value of 
E™ 1SS and/or of the scalar sum of reconstructed objects: leptons, jets and/or E™ 1SS is needed 


54 


For the chosen signal many neutrinos (> 6 ) appearing in the final state from X 4 and from r decay 
give rise to moderately high —when compared to signals from RPC supersymmetry— missing 
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Figure 5. Elf lss distribution for an LHC 
energy of y/s = 8 TeV and an integrated 
luminosity C = 20 fb -1 for a selected 
fi z/SSM point [53] . 


Figure 6. H\ versus lepton multiplicity for an 
LHC energy of y/s — 8 TeV and an integrated 
luminosity C — 20 fb -1 for a selected /iz/SSM 
point [53] . 


transverse energy, Elf lss , as depicted in Fig. [ 5 } Besides E™ 1SS , the scalar sum of the px of all 
reconstructed leptons, is also large in such events, as shown in Fig. [6j Alternatively, the 
sum of E™ 1SS and H ^ can be deployed for further background rejection. These observables can 
provide additional handles when selecting events with many leptons. In addition, the invariant 
masses, m^-i- and m 2 £+ 2t~ ma y prove useful for the purpose of signal distinction. 



Decay length (m) 


Figure 7. X 4 decay-length distribution 
for an LHC energy of y/s — 8 TeV and an 
integrated luminosity C — 20 fb -1 for a 
selected /xz/SSM point [53]. 



Figure 8. Charged-track mass versus the 
number of charge particles in each vertex for an 
LHC energy of y/s — 8 TeV and an integrated 
luminosity C — 20 fb -1 for a selected /iz/SSM 
point [53 . 


In the benchmark scenario under study, the X 4 is characterised by a proper lifetime of the 
order of r^o ~ 10 -9 s, which corresponds to a proper decay length of cr ~0 ~ 30 cm. This 
feature —quantified in Fig. [7| were the decay-length distribution is plotted— gives rise to 
displaced vertices. In a significant percentage of events, the y 9 would decay inside the inner 
tracker of the LHC experiments, e.g. in 28% of events it would decay within 30 cm and in 44% 
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events within 1 m. Therefore, the fiv SSM signal events would be characterised by displaced 
r leptons plus neutrinos. This distinctive signature opens up the possibility to exploit current 
or future variations of analyses carried out by ATLAS and CMS looking for a displaced muon 
and tracks 127,28] or searching for displaced dileptons [28][55], dijets [561 or muon jets [57 arising 
in Higgs decays to pairs of long-lived invisible particles. 

The kinematics of the DVs and their products in the chosen /iv SSM benchmark point have 
been studied thoroughly in Ref. |53|. The spacial distribution of DVs shows that an appreciable 
fraction of them falls in the inner-tracker volume of the LHC experiments, i.e. pov < 1 m and 
|^dv| < 2.5 m, thus DVs arising in the /ivSSM should be detectable at LHC, in principle, either 
with existing analyses [27, 28, 55-571 or via variations of those looking for displaced taus and 


£“ 1SS . The average x\ boost, on the other hand, as expressed by /3y, where /3 is X 4 velocity 
over c and 7 the Lorentz factor, is comparable to the signal analysed in an ATLAS related 
search for a muon and tracks originating from DVs [27 . The boost affects the efficiency with 
which such a DV can be reconstructed, since high y^boost leads to collimated tracks difficult 
to differentiate from primary vertices. In Fig. [ 8 j the correlation between the number of charged 
tracks in each DV, 3V tr k, and their invariant mass, mov, is shown. It has been demonstrated |27| 
that a selection of high-3V tr k and high-mov efficiently suppresses background from long-lived SM 
particles (B mesons, kaons). The modulation observed in N tr ^ is due to the one-prong or three- 
prong hadronic r decays. 


3.2. Higgs decays 

As pointed out earlier, /ivSSM can easily predict a Higgs boson of a mass compatible with the 
observed one 23, 24 . As becomes evident from Fig. [ 9 ] a 125-GeV Higgs boson, S 4 , can be 
accommodated within a wide range of tan/3 and A values |58l. 



tan (3 


Figure 9. Variation of m^ ree upper 
bound with tan/3 for different A values in 
the fivSSM. The horizontal lines represent 
the experimentally measured masses of the 
Higgs, m/i, and Z boson, mz- In the 
inset the region of tan/3 < 5 is shown for 


A <0.7 58 



Measured value 

m s o (GeV) 

ATLAS 


°- 2 -a6 ■ 

125.5 

/V+ 7 -CS 4 ) 

i-4i£5 ■ 

125.5 

M-n(SS) 

1.17 ± 0.27H 

125.4 

^ww*(S 4 ) 

. n +U.32 _ 

1 ' U -0.29 ■ 

125.5 

l‘zz-(S { !) 

1-441^ ■ 

125.36 

Combined 

l-30^;}« ■ 

125.5 

CMS 

m (S2) 

0.93 ± 0.49 ■ 

125 

Hr+r-(S° 4 ) 

0.91 ±0.27 ■ 

125 


1.13 ± 0.24 ■ 

125 


0.83 ± 0.21 ■ 

125 

Hzz-(S 4 ) 

1.00 ± 0.29 ■ 

125 

Combined 

1.00 ± 0.13 ■ 

125 


| ATLAS-CONF-2014-009 | 1408.7084 | 1408.5191 ■ CMS-PAS-HIG-14-009 


Figure 10. Measured signal strengths with 
la uncertainties obtained from the ATLAS and 
CMS experiments together with the correspond¬ 
ing values of Higgs mass [58]. 


This scalar would produce new states through two-body decays as long as these states are 
kinematically allowed. Decays to SfSj, P/’Pj 3 , X^ + 3 Xj +3 with i, j = 1,2,3 have been studied 
in detail in Ref. [58]. The question is whether these new decay widths respect the current 
measurements of the Higgs decays signal strength fixx- It is shown that the final states are 
dominated by a combination of prompt or displaced leptons, taus, jets and photons plus E^f lss 
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due to neutrinos. When the expected productions rates for these new signatures are compared 
with the current fixx measurements listed in Fig.[R)J it is found that they are compatible with 
the measured Higgs signal strengths fixx under the following conditions.: 

0.01 < A < 0 . 1 : All t^xx remain within 2 a of CMS measurements for 2.5 < tan/3 < 3.9. 

0.1 < A < 0.7 : Only (invisible) S 4 decays to x?+ 3 Xj +3 remain viable in the whole range of A. 

A > 0.1 : For decays to pair of binos, all /xxx are within 2a for 2.4 < tan/3 < 3.8. 

Hence there is plenty of room for new Higgs decays within the context of the “/x from z/’ 
supersymmetric standard model. 

3.3. New decays to Z/W 

In addition to the new SM-Higgs-like decays, /xz/SSM also introduces novel on-shell decays of the 
Z and W ± gauge bosons [591. These modes are typically encountered in regions of the parameter 
space populated with light singlet-like scalars, pseudoscalars and neutralinos. The complete 
spectrum of possible final states and their origin is presented in Table [T} The delayed “objects” 
occur in delayed decays of the neutralino, whereas the (short-lived) scalars and pseudoscalars 
deliver prompt products. 


Table 1 . Final states from non-standard Z and W ± decays with their respective origins 59 
The notation applied is x : e, /x, r, 7 , q and P, D stand for prompt and delayed, respectively. 


Z decay 

W ± decay 

2x D 2x D + E;f lss (via Xi+3Xj+3 ) 
2x p 2x p (via SfPj) 

£ p /t p + x D x D + Ef ss (via xfxj+3) 


These rare new decays are strongly constrained by the measurements of the Z and W total 
widths. Signatures with r-leptons and/or 6 -jets and with displaced objects would be preferred 
when probing these decays in order to suppress the huge SM background. Their low production 
rate, e.g. BR ~ 0(1O -5 ) for the Z, necessitates the high statistics that will become available 
when upcoming collides such as the GigaZ and TeraZ modes of the Linear Collider and TLEP, 
respectively, will be in operation |59|. 


4. Implications for dark matter 

We address here the issue of (not necessarily cold) dark matter in SUSY models with R -parity 
violation. It has been shown that these seemingly incompatible concepts can be reconciled in 
bRPV models with a gravitino [5j||| or an axino [l3[ LSP with a lifetime exceeding the age of 
the Universe. In both cases, RPV is induced by bilinear terms in the superpotential that can 
also explain current measurements on neutrino masses and mixings without invoking any GUT- 
scale physics. Decays of the next-to-lightest superparticle occur rapidly via RPV interactions, 
thus they do not disturb the Big-Bang nucleosynthesis, unlike the R -parity conserving case 60 


Decays of the NLSP into the gravitino and Standard Model particles do not contribute to 
the gravitino relic density in scenarios with broken R parity. This is because decay processes 
involving a gravitino in the final state are suppressed compared to R -parity breaking decays 
unless the amount of R -parity violation is extremely small. 

Interestingly, gravitino decays may produce monochromatic photons via G —» yz/, opening 
therefore the possibility to test this scenario with astrophysical searches, as we will see below. 
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Requiring the model parameters to correctly account for observed neutrino oscillation parameters 
implies that expected rates for 7 -ray lines produced by gravitino decays of mass above a few GeV 
would be confronted with the Fermi-LAT satellite observations, leading to an upper bound on 
the gravitino DM mass. For instance, the bRPV parameter region allowed by 7 -ray line searches, 
dark matter relic abundance, and neutrino oscillation data, has been determined obtaining a limit 
on the gravitino mass tUq <1 — 10 GeV corresponding to a relatively low reheat temperature |71. 
The allowed region in the gravitino mass versus gravitino lifetime plane is shown in Fig. m The 
yellow region is excluded by 7 -ray line searches: Fermi and EGRET constraints are above and 
below 1 GeV, respectively. 



Figure 11. Allowed gravitino mass- 
versus-lifetime region (grey color) con¬ 
sistent with neutrino oscillation data and 
astrophysical bounds on 7 -ray lines from 
dark matter decay for a bilinear RPV 
model. The lower and upper black 
lines correspond to mp 2 = 240 and 
3000 GeV, respectively. Adapted from 
Ref. [7]. 



Figure 12. Constraints on lifetime versus 
mass for decaying gravitino DM in the fiisSSM. 
Diagonal grey lines (shaded band) denote allowed 
(favoured) photino-neutrino mixing parameters for 
the fi z/SSM, The blue shaded area is excluded 
by the limits set by 5-year Fermi-LAT data. 
The upper bounds set by several other 7 -ray 
observations are also overlaid. From |5T|. 


Evidence on the four-year Fermi data that have found excess of a 130 GeV gamma-ray line 
from the Galactic Centre (GC) 1 61 1 have been studied in the framework of R- parity breaking 
SUSY. A decaying axino DM scenario based on the SUSY KSVZ axion model with the bilinear im¬ 
parity violation explains the Fermi 130-GeV gamma-ray line excess from the GC while satisfying 
other cosmological constraints |l4j. On the other hand, gravitino dark matter with trilinear 
RPV —in particular models with the LLE RPV coupling— can account for the gamma-ray 
line, since there is no overproduction of anti-proton flux, while being consistent with big-bang 
nucleosynthesis and thermal leptogenesis 1111 . 

Measurements of the cosmic-ray antiproton flux by the PAMELA experiment have been 
used recently to constrain the (decaying) gravitino mass and lifetime in the channels Zv, W£ 
and his |9[. Subsequently upper limits have been set on the size of the R-parity violating 
coupling in the range of 10 8 to 8 x 10 13 . These limits turn out to be more stringent than 
those coming from contributions to neutrino masses or from the baryon asymmetry in the early 
Universe. Combining them with lower limits from big-bang nucleosynthesis constraints on the 
NLSP lifetime allows narrowing down the available parameter space for gravitino DM with 
bRPV. 

Recent analyses of combined spectra of multiple galaxy clusters and the Andromeda galaxy 
from the XMM-Newton telescope, have revealed a tentative line with the central energy of 
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3.5 keV 1 62 1. First and foremost, any long lived particle that produces enough number of 
photons should be a good candidate as a source of X-rays. Among the various theoretical 
scenarios attempting to explain it, RPV SUSY with an LSP decaying to photon and a neutrino 
is a possibility. In bilinear RPV it was found that a warm-dark-matter axino with a mass of 
iria — 7 keV can have the proper lifetime and relic density to account for the observed X-ray 
emission line through its decay a —)> ju |63|. Alternatively the line may be due to annihilating 
DM with mu — 3.5 keV. Apart from the axino, other possible candidates include a gravitino, a 
bino, or a hidden sector photino as decaying DM leading to such a signal studied in the context 
of other-than-bilinear RPV SUSY scenarios [64] . 

Such gravitino DM is also proposed in the context of /ivSSM with profound prospects for 
detecting 7 rays from their decay |49 . In studies on the prospects of the Fermi-LAT telescope 
to detect such monochromatic lines in five years of observations of the most massive nearby 
extragalactic objects, it was found that a gravitino with a mass range of 0.6 — 2 GeV, and 
with a lifetime range of about 3 x 10 27 — 2 x 10 28 s should be detectable with a signal-to-noise 
ratio larger than three [501. After confronting the actual 7 -ray flux data, limits on the model 
parameters have been set|50,511. It was shown that gravitino masses larger than about 4 GeV 
are already excluded in the fin SSM by Fermi-LAT data of the galactic halo 1501. The constraints 
on G lifetime versus mass obtained by the 5-year Fermi-LAT and other 7 -ray observations, as a 
consequence of the limits on line emission, are summarised in Fig. [12] If the photino-neutrino 
mixing parameter \U^ V \ from the 7x7 neutralino mixing matrix is between 10 -16 and 10 —12 , 
the correct neutrino masses are reproduced in /iz/SSM. Values of the gravitino mass larger than 
5 GeV, are disfavoured, or lifetimes smaller than ~ 10 28 s, are excluded at 95% CL as DM 
candidates. 

R -parity breaking couplings can be sufficiently large to lead to interesting expectations for 
collider phenomenology. The neutralino NLSP, depending on the RPV model considered, may 
decay into 6.65] 

Xi Xi G 7 , 

Xi -> Xi -> GZ, (7) 

X? v li Xi Gh. 

Such decays may be probed at the LHC via inclusive channels characterised by leptons, 

many jets, large E™ 1SS and/or photons or exclusive channels involving the reconstruction of a 
Z or a h from its decay products. The cases where the NLSP is long-lived yet with a decay 
length comparable to the dimensions of an LHC experiment are particularly interesting, as they 
give rise to displaced tracks/leptons and non-pointing photons. The possibility to measure the 
neutralino decay length provides an extra handle to constrain the underlying SUSY model. 
Nevertheless determining whether R parity is conserved or broken may not be trivial since the 
DM particles themselves, whether absolutely stable or quasi-stable, cannot be detected directly 
in collider experiments. To this effect, an interplay between collider and astroparticle physics is 
necessary to pin down the dark matter properties and the related theoretical-model parameters. 

5. Summary and outlook 

The hitherto null results of searches for supersymmetry in conventional channels calls for a 
more systematic and thorough consideration of non-standard SUSY theoretical scenarios and 
experimental techniques. To this effect, scenarios involving violations of R parity and/or 
(meta)stable particles arise as interesting alternatives. 

R-parity violating supersymmetry may reproduce correctly the measured neutrino physics 
parameters. Moreover its enriched mass spectrum and R-parity breaking decays can lead to 
novel signals at colliders, among which very few possibilities have been explored so far. In 



12 


V A Mitsou 


particular, searches for displaced objects at the LHC offer an attractive possibility as well as 
the prompt multilepton analyses. The former suffer from less background sources from SM 
processes when compared to prompt-object searches. More sophisticated searches are expected 
by ATLAS, CMS and LHCb in the near future during LHC Run II, where energies of 13 — 14 TeV 
will be probed. 

Acknowledgments 

The author is grateful to the DISCRETE 2014 organising committee for the kind invitation 
to the Symposium. She thanks Pradipta Ghosh for providing material for the paper. She 
acknowledges support by the Spanish Ministry of Economy and Competitiveness (MINECO) 
under the project FPA2012-39055-C02-01, by the Generalitat Valenciana through the project 
PROMETEO 11/2013-017 and by the Spanish National Research Council (CSIC) under the 
JAE-Doc program co-funded by the European Social Fund (ESF). 

References 

[1] Nilles H P 1984 Phys. Rept. 110 1; 

Martin S P 1998 A Supersymmetry primer (Preprint hep-ph/9709356) In Perspectives on supersymmetry II 
ed G L Kane (Singapore: World Scientific) 

[2] For a review on dark-matter searches in colliders, see e.g.: Mitsou V A 2013 Int. J. Mod. Phys. A 28 1330052 

(Preprint arXiv: 1310.1072 [hep-ex]) 

[3] Evans L and Bryant P 2008 J. Instrum. JINST03(2008)S08001 

[4] Barbier R et al 2005 Phys. Rept. 420 1 (Preprint hep-ph/0406039) 

[5] Takayama F and Yamaguchi M 2000 Phys. Lett. B 485 388 (Preprint hep-ph/0005214) 

[6] Hirsch M, Porod W and Restrepo D 2005 J. High Energy Phys. JHEP03(2005)062 (Preprint hep-ph/0503059) 

[7] Restrepo D, Taoso M, Valle J W F and Zapata O 2012 Phys. Rev. D 85 023523 (Preprint arXiv:1109.0512 

[hep-ph]) 

[8] Cottin G, Diaz M A, Guzman M J and Panes B 2014 Eur. Phys. J. C 74 3138 (Preprint arXiv: 1406.2368 

[hep-ph]) 

[9] Delahaye T and Grefe M 2013 J. Cosmol. Astropart. Phys. JCAP12(2013)045 (Preprint arXiv:1305.7183 

[hep-ph]) 

[10] Steffen F D 2009 Eur. Phys. J. C 59 557 (Preprint arXiv:0811.3347 [hep-ph]) 

[11] Liew S P 2013 Phys. Lett. B 724 88 (Preprint arXiv: 1304.1992 [hep-ph]) 

[12] Holman R, Lazarides G and Shah Q 1983 Phys. Rev. D 27 995; 

Nomura Y and Thaler J 2009 Phys. Rev. D 79 075008 (Preprint arXiv:0810.5397 [hep-ph]); 

Baer H, Kraml S, Lessa A and Sekmen S 2010 J. Cosmol. Astropart. Phys. JCAP11(2010)040 (Preprint 
arXiv: 1009.2959 [hep-ph]); 

Baker K et al 2013 Ann. Phys. 525 A93 (Preprint arXiv: 1306.2841 [hep-ph]); 

Bae K J, Baer H and Chun E J 2013 J. Cosmol. Astropart. Phys. JCAP12(2013)028 (Preprint 
arXiv: 1309.5365 [hep-ph]); 

Graf P and Steffen F D J. Cosmol. Astropart. Phys. JCAP12(2013)047 (Preprint arXiv: 1302.2143 [hep-ph]) 

[13] Chun E J and Kim H B 2006 J. High Energy Phys. JHEP10(2006)082 (Preprint hep-ph/0607076) 

[14] Endo M, Hamaguchi K, Liew S P, Mukaida K and Nakayama K 2013 Phys. Lett. B 721 111 (Preprint 

arXiv: 1301.7536 [hep-ph]) 

[15] Porod W, Hirsch M, Romao J and Valle J W F 2001 Phys. Rev. D 63 115004 (Preprint arXiv:hep-ph/0011248) 

[16] Romao J C, Diaz M A, Hirsch M, Porod W and Valle J W F 2000 Phys. Rev. D 61 071703 (Preprint 

arXiv:hep-ph/9907499) 

[17] Hirsch M, Diaz M A, Porod W, Romao J C and Valle J W F 2000 Phys. Rev. D 62 113008 [Erratum-ibid. 

2002 D 65 119901] (Preprint arXiv:hep-ph/0004115) 

[18] Diaz M A, Hirsch M, Porod W, Romao J C and Valle J W F 2003 Phys. Rev. D 68 013009 [Erratum-ibid. 

2005 D 71 059904] (Preprint arXiv:hep-ph/0302021) 

[19] Abada A and Losada M 2000 Phys. Lett. B 492 310 (Preprint arXiv:hep-ph/0007041) 

[20] de Campos F, Eboli O J P, Magro M B, Porod W, Restrepo D, Das S P, Hirsch M and Valle J W F 2012 

Phys. Rev. D 86 075001 (Preprint arXiv: 1206.3605 [hep-ph]) 

[21] de Campos F, Eboli O J P, Magro M B, Porod W, Restrepo D, Hirsch M and Valle J W F 2008 J. High 

Energy Phys. JHEP05(2008)048 (Preprint arXiv:hep-ph/0712.2156) 


Exploring neutrino physics at LHC via R-parity violating SUSY 


13 


[22] Porod W 2003 Comput. Phys. Commun. 153 275 (Preprint hep-ph/0301101); 

Porod W and Staub F 2012 Comput. Phys. Commun. 183 2458 (Preprint arXiv: 1104.1573 [hep-ph]) 

[23] Aad G et al [ATLAS Collaboration] 2012 Phys. Lett. B 716 1 ( Preprint arXiv: 1207.7214 [hep-ex]) 

[24] Chatrchyan S et al [CMS Collaboration] 2012 Phys. Lett. B 716 30 (Preprint arXiv: 1207.7235 [hep-ex]) 

[25] Torro E, Mitsou V A and Garcia C 2010 Prospects for detecting imparity violating supersymmetry with 

ATLAS in the muon plus jets channel ATLAS Note CERN-ATL-PHYS-INT-2010-042 

[26] Torro Pastor E 2013 Study of supersymmetric signals with 17-parity violation in ATLAS at LHC PhD thesis 

CERN-THESIS-2012-308 (Univ. Valencia) 

[27] Aad G et al [ATLAS Collaboration] 2013 Phys. Lett. B 719 280 (Preprint arXiv: 1210.7451 [hep-ex]); 

Aad G et al [ATLAS Collaboration] 2012 Phys. Lett. B 707 478 (Preprint arXiv: 1109.2242 [hep-ex]) 

[28] Aad G et al [ATLAS Collaboration] 2015 Search for massive, long-lived particles using multitrack displaced 

vertices or displaced lepton pairs in pp collisions at y/s — 8 TeV with the ATLAS detector at the LHC 
Preprint arXiv: 1504.05162 [hep-ex] submitted to Phys. Rev. D 

[29] de Campos F, Eboli O J P, Magro M B and Restrepo D 2009 Phys. Rev. D 79 055008 (Preprint 

arXiv:0809.0007 [hep-ph]); 

de Campos F, Eboli O J P, Magro M B, Restrepo D and Valle J W F 2009 Phys. Rev. D 80 015002 
(Preprint arXiv:0809.1637 [hep-ph]); 

de Campos F, Eboli O J P, Magro M B, Porod W, Restrepo D and Valle J W F 2005 Phys. Rev. D 71 
075001 (Preprint hep-ph/0501153) 

[30] De Campos F, Diaz M A, Eboli O J P, Magro M B and Mercadante P G 2002 Nucl. Phys. B 623 47 (Preprint 

hep-ph/0110049); 

de Campos F, Diaz M A, Eboli O J P, Lineros R A, Magro M B and Mercadante P G 2005 Phys. Rev. D 
71 055008 (Preprint hep-ph/0409043); 

De Campos F, Diaz M A, Eboli O J P, Magro M B, Porod W and Skadhauge S 2008 Phys. Rev. D 77 
115025 (Preprint arXiv:0803.4405 [hep-ph]) 

[31] Aad G et al [ATLAS Collaboration] 2012 Phys. Rev. D 85 012006 (Preprint arXiv: 1109.6606 [hep-ex]) 

[32] ATLAS Collaboration 2012 Search for supersymmetry at y/s m 7 TeV in final states with large jet multiplicity, 

missing transverse momentum and one isolated lepton with the ATLAS detector ATLAS Note 
ATLAS-CONF-2012-140 http://cds.cern.ch/record/1483511 

[33] Mitsou V A 2014 EPJ Web Conf. 71 00092 (Preprint arXiv: 1401.2573 [hep-ph]) 

[34] Aad G et al [ATLAS Collaboration] 2014 J. High Energy Phys. JHEP06(2014)035 (Preprint arXiv: 1404.2500 

[hep-ex]) 

[35] Aad G et al [ATLAS Collaboration] 2014 J. High Energy Phys. JHEP09(2014)103 (Preprint arXiv: 1407.0603 

[hep-ex]) 

[36] Aad G et al [ATLAS Collaboration] 2015 J. High Energy Phys. JHEP04(2015)116 (Preprint arXiv: 1501.03555 

[hep-ex]) 

[37] Khachatryan V et al [CMS Collaboration] 2015 J. High Energy Phys. JHEP04(2015)124 (Preprint 

arXiv: 1502.06031 [hep-ex]) 

[38] Aad G et al [ATLAS Collaboration] 2015 Search for supersymmetry in events containing a same-flavour 

opposite-sign dilepton pair, jets, and large missing transverse momentum in y/s = 8 TeV pp collisions with 
the ATLAS detector Preprint arXiv: 1503.03290 [hep-ex] Eur. Phys. J. C to appear 

[39] Papucci M, Ruderman J T and Weiler A 2012 J. High Energy Phys. JHEP09(2012)035 (Preprint 

arXiv: 1110.6926 [hep-ph]) 

[40] Djouadi A et al [MSSM Working Group Collaboration] 1999 The Minimal supersymmetric standard model: 

Group summary report Preprint hep-ph/9901246 

[41] Cahill-Rowley M W, Hewett J L, Ismail A and Rizzo T G 2012 Phys. Rev. D 86 075015 (Preprint 

arXiv: 1206.5800 [hep-ph]) 

[42] ATLAS Collaboration 2015 Constraints on promptly decaying supersymmetric particles with lepton-number- 

and P-parity-violating interactions using Run-1 ATLAS data ATLAS Note ATLAS-CONF-2015-018 
http://cds.cern.ch/record/2017303 

[43] Lopez-Fogliani D E and Munoz C 2006 Phys. Rev. Lett. 97 041801 (Preprint hep-ph/0508297) 

[44] Escudero N, Lopez-Fogliani D E, Munoz C and Ruiz de Austri R 2008 J. High Energy Phys. JHEP12(2008)099 

(Preprint arXiv:0810.1507 [hep-ph]) 

[45] Kim J E and Nilles H P 1984 Phys. Lett. B 138 150 

[46] Ellwanger U, Hugonie C and Teixeira A M 2010 Phys. Rept. 496 1 (Preprint arXiv:0910.1785 [hep-ph]) 

[47] Nilles H-P and Polonsky N 1997 Nucl. Phys. B 484 33 (Preprint hep-ph/9606388) 

[48] Ghosh P and Roy S 2009 J. High Energy Phys. JHEP04(2009)069 (Preprint arXiv:0812.0084 [hep-ph]); 

Fidalgo J, Lopez-Fogliani D E, Munoz C and Ruiz de Austri R J. High Energy Phys. JHEP10(2011)020 
(Preprint arXiv: 1107.4614 [hep-ph]); 


14 


V A Mitsou 


Ghosh P, Dey P, Mukhopadhyaya B and Roy S 2010 J. High Energy Phys. JHEP05(2010)087 ( Preprint 
arXiv: 1002.2705 [hep-ph]); 

Bartl A, Hirsch M, Vicente A, Liebler S and Porod W 2009 J. High Energy Phys. JHEP05(2009)120 
{Preprint arXiv:0903.3596 [hep-ph]); 

Bandyopadhyay P, Ghosh P and Roy S 2011 Phys. Rev. D 84 115022 {Preprint arXiv:1012.5762 [hep-ph]); 
Liebler S and Porod W 2012 Nucl. Phys. B 855 774 {Preprint arXiv: 1106.2921 [hep-ph]); 

Fidalgo J, Lopez-Fogliani D E, Munoz C and Ruiz de Austri R 2009 J. High Energy Phys. JHEP08(2009)105 
{Preprint arXiv:0904.3112 [hep-ph]); 

Fidalgo J and Munoz C 2012 J. High Energy Phys. JHEP04(2012)090 {Preprint arXivill 11.2836 [hep-ph]) 

[49] Choi K-Y, Lopez-Fogliani D E, Munoz C and Ruiz de Austri R 2010 J. Cosmol. Astropart. Phys. 

JCAP03(2010)028 {Preprint arXiv:0906.3681 [hep-ph]) 

[50] Gomez-Vargas G A, Fornasa M, Zandanel F, Cuesta A J, Munoz C, Prada F and Yepes G 2012 J. Cosmol. 

Astropart. Phys. JCAP02(2012)001 {Preprint arXiv: 1110.3305 [astro-ph.HE]) 

[51] Albert A et al [Fermi-LAT Collaboration] 2014 J. Cosmol. Astropart. Phys. JCAP10(2014)023 {Preprint 

arXiv: 1406.3430 [astro-ph.HE]) 

[52] Chung D J H and Long A J 2010 Phys. Rev. D 81 123531 {Preprint arXiv: 1004.0942 [hep-ph]) 

[53] Ghosh P, Lopez-Fogliani D E, Mitsou V A, Munoz C and Ruiz de Austri R 2013 Phys. Rev. D 88 015009 

{Preprint arXiv: 1211.3177 [hep-ph]) 

[54] Chatrchyan S et al [CMS Collaboration] 2014 Phys. Rev. D 90 032006 {Preprint arXiv: 1404.5801 [hep-ex]); 

Aad G et al [ATLAS Collaboration] 2014 Phys. Rev. D 90 052001 {Preprint arXiv: 1405.5086 [hep-ex]) 

[55] Khachatryan V et al [CMS Collaboration] 2015 Phys. Rev. D 91 052012 {Preprint arXiv: 1411.6977 [hep-ex]) 

[56] Khachatryan V et al [CMS Collaboration] 2015 Phys. Rev. D 91 012007 (Preprint arXiv: 1411.6530 [hep-ex]); 

Aad G et al [ATLAS Collaboration] 2015 Phys. Lett. B 743 15 {Preprint arXiv: 1501.04020 [hep-ex]) 

[57] Aad G et al [ATLAS Collaboration] 2013 Phys. Lett. B 721 32 {Preprint arXiv: 1210.0435 [hep-ex]) 

[58] Ghosh P, Lopez-Fogliani D E, Mitsou V A, Munoz C and Ruiz de Austri R 2014 J. High Energy Phys. 

JHEP11(2014) 102 {Preprint arXiv: 1410.2070 [hep-ph]) 

[59] Ghosh P, Lopez-Fogliani D E, Mitsou V A, Munoz C and Ruiz de Austri R 2015 Phys. Rev. D 91 035020 

{Preprint arXiv: 1403.3675 [hep-ph]) 

[60] Fukugita M and Yanagida T 1986 Phys. Lett. B 174 45; 

Buchmuller W, Covi L, Hamaguchi K, Ibarra A and Yanagida T 2007 J. High Energy Phys. 
JHEP03(2007)037 {Preprint hep-ph/0702184); 

Cyburt R H, Ellis J, Fields B D, Luo F, Olive K A and Spanos V C 2009 J. Cosmol. Astropart. Phys. 
JCAP10(2009)021 {Preprint arXiv:0907.5003 [astro-ph.CO]) 

[61] Bringmann T, Huang X, Ibarra A, Vogl S and Weniger C 2012 J. Cosmol. Astropart. Phys. JCAP07(2012)054 

{Preprint arXiv: 1203.1312 [hep-ph]); 

Weniger C 2012 J. Cosmol. Astropart. Phys. JCAP08(2012)007 {Preprint arXiv: 1204.2797 [hep-ph]); 
Tempel E, Hektor A and Raidal M 2012 J. Cosmol. Astropart. Phys. JCAP09(2012)032 [Addendum-ibid. 
JCAP11(2012)A01] {Preprint arXiv: 1205.1045 [hep-ph]) 

[62] Bulbul E, Markevitch M, Foster A, Smith R K, Loewenstein M and Randall S W 2014 Astrophys. J. 789 13 

{Preprint arXiv: 1402.2301 [astro-ph.CO]); 

Boyarsky A, Ruchayskiy O, Iakubovskyi D and Franse J 2014 Phys. Rev. Lett. 113 251301 {Preprint 
arXiv:1402.4119 [astro-ph.CO]) 

[63] Park J C, Park S C and Kong K 2014 Phys. Lett. B 733 217 {Preprint arXiv: 1403.1536 [hep-ph]); 

Choi K Y and Seto O 2014 Phys. Lett. B 735 92 {Preprint arXiv: 1403.1782 [hep-ph]) 

[64] Kolda C and Unwin J 2014 Phys. Rev. D 90 023535 {Preprint arXiv: 1403.5580 [hep-ph]); 

Bomark N-E and Roszkowski L 2014 Phys. Rev. D 90 011701 {Preprint arXiv: 1403.6503 [hep-ph]); 

Liew S P 2014 J. Cosmol. Astropart. Phys. JCAP05(2014)044 (Preprint arXiv: 1403.6621 [hep-ph]) 

[65] Hajer J 2013 Long-lived neutralinos as probes of gravitino dark matter PhD thesis DESY-THESIS-2013-031 

(Univ. Hamburg); 

Bobrovskyi S, Hajer J and Rydbeck S 2013 J. High Energy Phys. JHEP02(2013)133 {Preprint 
arXiv: 1211.5584 [hep-ph]) 


